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1 Introduction and summary 



In the theory of phase transitions, the relevance of the analytic continuation of the free energy to 
complex values of its parameters is known since the work of Lee and Yang [1]. They have shown that 
the thermodynamical equation of state is completely determined by the distribution of zeros of the 
partition function in the fugacity complex plane, or in the magnetic field complex plane in the case of 
spin systems. The distribution of zeros determine the analytical properties of the free energy - they 
typically accumulate in arrays and produce branch cuts. From the analytical properties of the free 
energy one can reconstruct the qualitative picture of the critical phases and the phase boundaries. 

The most studied example is the ferromagnetic Ising model in a complex magnetic field. In 
the high-temperature phase (T > Tent), the free energy as a function of the magnetic field has two 
branch cuts, starting at the points = ±.iHc{T) on the imaginary axis [2]. The gap between the 
two branch points vanishes at the critical temperature Tent. Fisher [3] named the branch points Lee- 
Yang singularities, and suggested that they are described by a non-unitary field theory with a cubic 
potential. Cardy [4] identified the Lee- Yang singularity as the simplest non-unitary CFT, the minimal 
model M.2,b with a central charge = —22/5. On the other hand, in the low-temperature phase 
(T < Tcrit), the free energy has a weak singularity at the origin, predicted by Langer's theory [5,6]. 
The analytical continuation from the positive part of the real axis has a branch cut on the negative axis, 
known as Langer branch cut. The discontinuity across this cut is exponentially small when T — > Tcrit. 
It can be explained with condensation of droplets of the stable phase near a metastable point. 

A thorough analysis of the analytic properties of the free energy near the critical point {H = 
0, T = Tcrit), based on the truncated CFT approach [7], was performed by Fonseca and Zamolod- 
chikov [8]. In the continuum limit, the Ising model is described by an euclidean quantum field theory, 
the Ising field theory (IFT), which is a double perturbation of massless Majorana fermion. The Ising 
field theory is formally defined by the action 

•AiFi = Ac=i/2 ~ t Js{x)cfx — h Js{x)d^x, (1.1) 

where the first term is the action of a free massless Majorana fermion, s{x) is the spin field with 
conformal weights As = As = 1/16 and e is the energy density with conformal weight A^ = 
Ag = 1/2. The coupling constants t and h are the renormalized temperature and magnetic field, they 
parametrize the vicinity of the critical point T = Tent, H = 0. In the QFT language, the specific free 
energy is interpreted as the vacuum energy density, and its discontinuity along the Langer branch cut 
as the decay rate of a "false vacuum" [9, 10]. By dimensional arguments, the vacuum energy density 
of the Ising field theory must have the form 

F{t,h)^t''logt'' + t^G{0, (1.2) 

where the first term is the free energy of a (massive) Majorana fermion and G{^) is a scaling function 
of the dimensionless strength of the magnetic field, 

? = (1-3) 

The high temperature (t > 0) and the low-temperature (t < 0) regimes are described by two different 
scaling functions, Ghigh(0 GiowiO- 

The function Ghigh(0 is analytic at small ^ and has two symmetric Yang-Lee branch points at 
^ = ±i^c on the imaginary axis (Fig. 1), while the function Giow(0 characterized by Langer 
branch cut starting at the origin (Fig. 2). The two scaling functions are analytically related in the 
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Figure 1: Analytic structure of G'high(C) 
in the complex ^-plane. The function has 
two singularities (YL edges) at ^ = ±iCc 
and a third one at£^ — oo (the critical point). 
There are two cuts along the imaginary axis 
connecting the two YL edges and the critical 
point. 



Figure 2: Analytic structure of Giow{£_). 
The function has two singularities, the criti- 
cal point at ^ = cx) and the low-temperature 
fixed point at ^ 0. The two singulari- 
ties are connected by a branch cut extending 
from — oo to 0. 



vicinity of the point ^ = oo because the magnetic field creates a mass gap and smears the phase 
transition at t = 0. Therefore Ghigh and Giow must correspond to two different branches of the same 
multivalued function of the variable ^. These analytic properties represent the standard analyticity 
assumptions; they do not exclude the presence of more singularities in the other sheets of G(^). 

The numerical results obtained in [8] led the authors to a stronger assumption, which they referred 
to as 'extended analyticity', and which states that the Yang-Lee edge is the nearest singularity un- 
der the Langer branch cut. The extended analyticity suggests an interpretation of the YL edge as a 
quantum counterpart of the spinodal point in the classical theory of phase transitions. The classical 
low-temperature free energy is regular at = 0, but shows a branch-cut singularity at some negative 
H = —Hsp, where the metastable phase becomes classically unstable. It was argued in [8] that the 
spinodal singularity does not disappear completely due to the quantum effects, but moves under the 
Langer cut and reappears as YL edge in the high-temperature phase. More recently, the numerical 
results of [8] were confirmed by a calculation based on the lattice formulation of the problem [1 1]. 

It is tempting to think that the qualitative picture of the analytic properties of the free energy, 
proposed in [8], applies to all statistical systems having a line of first order transitions ending at a 
second order transition point, and that the extended analyticity conjecture has a universal character. 
However, it is not clear how to address this issue directly. The Langer singularity is too weak to to be 
measured experimentally, and the TCFT approach used in [8] is sufficiently reliable only in the case 
of Majorana fermions. ^ 

On the other hand, it is known that some questions about the two-dimensional statistical systems 
can be answered by the following detour: formulate the system on a dynamical lattice, solve the 
problem exactly and then interpret the answer for the original system. This procedure is usually 
called 'coupling to 2D quantum gravity'. Due to the enhanced symmetry, which erases the coordinate 
dependence of the correlation functions, the systems on dynamical lattices are much simple to resolve. 
This approach is based on the fact that there is a well established correspondence between the universal 

'An attempt to use the approach of of [8] to study the tricritical Ising model in magnetic field was made in [12]. 
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properties on the flat and dynamical lattices. On the dynamical lattice, the description of the critical 
points is given in terms of Liouville gravity [13]. The matter in Liouville gravity is given by the CFT 
describing the critical point of the statistical system on a flat lattice, while the Liouville field describes 
the fluctuations of the metric. The matter fields influence the geometry in a precise way, which is taken 
into account by accompanying Liouville vertex operators [14-16]. Knowing the correlation functions 
in 2D gravity, one can reconstruct the scaling dimensions of the matter operators. Moreover, the 
method based on coupling to 2D gravity is efficient also away from the critical points, where it makes 
possible to reconstruct the bulk and the boundary renormalization flow diagrams of the matter theory 
on a flat lattice, knowing the exact solution on dynamical lattice. The corresponding field theory is 
Liouville gravity perturbed by one or several fields. 

The ising model is the only nontrivial example where the analytic properties of the free energy are 
known both on flat and dynamical lattices. The exact solution of Boulatov and Kazakov [17] shows 
that the Yang-Lee singularities appear also on a dynamical lattice and that near each YL singularity 
the critical behavior is that of a Yang-Lee CFT coupled to gravity [18-20]. However, it is not a 
priori clear if coupling to 2D gravity could be used to study the decay of the metastable vacuum by 
nucleation. Such a possibility was first explored in an unpublished work by Al. Zamolodchikov [21] 
and a subsequent work by A. and Al. Zamolodchikov [22]. 

The universal behavior of the Ising model on dynamical lattice is described by Ising Quantum 
Gravity (IQG), which represents a double perturbation of Liouville gravity with matter centtal charge 
1/2. The IQG is formally defined by the action 

^IQG = A=l/2 Liouville gravity " t J e''^^'*'^^^ <f X -h j a{x) ^'^'^^^^U^X. (1.4) 

The exponents of Liouville dressing operators, = 1/3 and = 5/6, are such that the dressed 
matter operators become densities [15, 16]. The vacuum energy density of IQG must be of the form 

:FiQG{t,h) = f'g{0. ^ = h/t'/\ (1.5) 

where the exponent 5/2 = a^/cKe is such that ^ is invariant under rescalings of the metric. The 
vacuum energy density of the IQG, or the universal specific free energy of the discrete model, was ex- 
tracted in [21] from the Bulatov-Kazakov exact solution [17]. For some normahzation of the couphng 
constants, the universal free energy is given in parametric form by 

:FiQG{t, h) = -\u''-\v?-lh''{u + , 

=u{u + t)\ ^ ^ ^ 

The scaling function Q{Cj defined in (1.5) is analytic on a Riemann surface representing a 5-sheet 
covering of the complex -plane. The Riemann surface has a high-temperature (HT) sheet with two 
YL branch points at ^ = ±i^c and two copies of the low-temperature (LT) sheet having a single cut 
starting at ^ = and going to infinity. The HT sheet is coimected to the two LT sheets by two auxiliary 
sheets. In this way, the scaling functions t/high and ^low are identified with the restrictions of a single 
function Q to the HT and the LT sheet of its Riemann surface. 

In IQG, the discontinuity of Q\o^ vanishes near the origin as ^^/^ and not exponentially, as pre- 
dicted by the droplet model [21]. A field-theoretical explanation of the power law was given in [22], 
where Langer's theory of the condensation point was adjusted to the case of a fluctuating metric. In the 
case of a flat lattice, the discontinuity across the Langer cut is related to the Boltzmann weight of the 
critical droplet of the true vacuum in the sea of the false vacuum (Fig.3). A. and Al. Zamolodchikov 
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argued in [22] that in presence of gravity the droplets of the stable phase influence the metric, so that 
the perimeter of the boundary of a droplet grows slower than the square of the area (Fig. 4). This 
leads to a much faster decay of the unstable vacuum, which changes the character of the singularity at 
^ = 0. This phenomenon was named by the authors of [22] 'critical swelling'. 

The work [22] gave a convincing evidence that the Langer singularity in IFT has its counterpart in 
IQG, and that both singularities are due to the presence of a weakly metastable state. Is it possible to 
extend the IFT/IQG correspondence to more general spin systems? In fact, the analysis of [22] applies 
to any droplet model in presence of gravity. Due to the critical swelling, the Boltzmann weight of the 
critical droplet is expected to behaves as a power, /i^'i"" , where the exponent giov, is determined by 
the conformal anomaly in the low-temperature phase Qow- The saddle point analysis of [22] gives 
5iow = — gCiow in the classical limit ciow — — oo, where the fluctuations of the Liouville field can 
be neglected. The authors of [22] conjectured (we will refer to this as Z-Z conjecture) that the exact 
nucleation exponent (710W is related to the central charge of the stable phase as follows, 

p 6 

5I0W = 7 for Clow = 1 7 7T- (1-7) 

p — 1 p{p — 1) 

In the case of IQG, ciow = and gio^ = 3/2, which is in agreement with the singularity of ^low at 
^ = 0. If the conjecture (1.7) is confirmed, this would open the possibility to study the analytical 
properties of the free energy of a spin systems on a flat lattice, in particular in what concerns the 
Langer and YL singularities, using its exact solution on a dynamical lattice. 

The conjecture (1.7) about the nucleation rate exponent can be tested by looking for more solvable 
examples among the dynamical lattice statistical systems, which exhibit metastability and have non- 
trivial massless modes in the stable vacuum. In this paper we argue that the 0{n) loop model in 
external magnetic field is such a system. 

Technically the main result of this paper is the explicit expression for the specific free energy 
of the gravitational 0(n) model in the continuum limit, which we obtain using the correspondence 
with a large N matrix model [23]. The main body of the text is devoted to the interpretation of 
this exact solution. In Sect. (2) we formulate the expected analytic properties of the free energy 
of the 0{n) loop model on a flat lattice as well as the the non-perturbative effects associated with 
the existence of a hypothetical metastable vacuum. In Sect. (3) we perform the same analysis of 
the gravitational 0(n) model. In particular, assuming the existence of a metastable state and the 
conjectured nucleation exponent (1.7), we reproduce the form of the series expansion of the free 
energy of the low-temperature phase at h = 0. In Sect. 4 we verify that the exact expression for the 
free energy possesses the expected analytic properties, which justifies the assumption of the existence 
of a weakly metastable state in the 0{n) loop model. Our result verify the Z-Z conjecture for the 



4 



Figure 5: A polymer configuration on a 
honeycomb lattice 



Figure 6: A polymer configuration on a pla- 
nar graph with the topology of a disk 



wliole spectrum — oo < c < 1 of the central charge of the matter field. Since the paper is somewhat 
technical, below we give a short summary of our results. 



o o o 



The 0(n) loop model [24, 25] has a geometrical formulation in terms of nonintersecting loops 
with fugacity n. In presence of magnetic field H, the geometrical expansion involves also open lines 
with fugacity H^, as in the example shown in Fig. 5. The temperature T controls the total length 
of the linear polymers. In the conformal window — 2 < n < 2, the theory has a conformal invariant 
critical point at T = Tc and H = 0, described by a CFT with central charge [26], 

c = 1 - r, (1-8) 

where p is related to the dimensionality n by 

n = 2cos-, p>l. (1.9) 
P 

The critical 0{n) model covers the whole one-parameter family of universality classes having mass- 
less modes in the low-temperature phase, with a central charge — oo < c < 1. 

The phase diagram of the 0{n) loop model is similar to that of the Ising model, which corresponds 
to the particular case n = 1, or p = 3. For T > Tc the typical linear polymer is short and the theory 
has a mass gap. Nienhuis [25] showed, using the mapping to other solvable lattice models and the 
Coulomb gas techniques, that for T < Tc and H = the model is characterized by a massless low- 
temperature (LT) phase, also known as a phase of dense loops. The LT phase of the 0{n) vector 
model is described by another CFT with central charge^ 

/? 

Clow = 1 - TV. (1.10) 

P{P - 1) 

The vicinity of the critical point (T = Tc, H = 0) is described by the 'O(n) field theory', or 
shortly 0„FT, which is a perturbation of the CFT with central chai^ge (1.8). The 0{n) field theory is 



"To avoid confusion, let us mention that this low-temperature phase is not a Goldstone phase. The 0(n) loop model 
flows to a Goldstone phase having central charge n — 1 appears when the loops are allowed to intersect with non-zero 
probability [27,28], and in the microscopic formulation given in [24,25] such intersections are strictly forbidden. 
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formally defined by the action 



Ao„FT = Arit - hj S{x) (fx-t J £{x) (fx , (1.11) 

where s and e are the conformal fields generating the perturbations respectively in the magnetic field 
h H and in the temperature t ^ T — T^. When p is an odd integer, p = 2m—l, this action describes 
the unitary minimal conformal theory M.2m-i,2m perturbed by the two primary fields 

£ = $1,3, S = ^m,m {m = \{p + I)) . (1.12) 

The Ising field theory corresponds to m = 2, or n = 1. The thermal flow generated by e ends, 
depending on the sign of the temperature coupUng, at a massive theory in the HT regime, or at a CFT 
with the lower central charge (1.10) in the LT regime [29]. For a general p > 1, the 0(n) field theory 
does not exist as a local field theory, but it has nonetheless an unambiguous definition in terms of 
grand canonical ensemble of linear polymers on the lattice. 

In Sect. 2 we speculate about the analytic properties of the specific free energy of OnFT, assuming 
that there exists a weakly metastable state and Langer singularity at least in some finite vicinity of 
the point p = 3. Since the free energy depends analytically on the parameter p, this is a natural 
assumption, although for p 7^ 3 we do not see any geomettically ttansparent picture of the nucleation 
mechanism.^ We give a heuristic argument about the form of the Langer singularity for p 7^ 3. The 
non-perturbative corrections should be of the form 

i^nonpert ~ / exp(-COnSt//), / ~ (1.13) 

where A'"^ is the conformal dimension of the 0{n) spin field s in the low-temperature phase. 

As in the case of IFT, the specific free energy is determined by two different seating functions in 
the HT and in the LT regimes:^ 



F{t,h) 



' t"^ Ghigh(e), i = t-^^h if t>0, 



(1.14) 



where the power 2±i and the exponent 

, . _ 1 - As (3p + l)(5p + 3) 

= 1^ = 32^ ^^-^^^ 

are determined by the conformal weights = A^ and Ag = Ag respectively of e and s at the critical 
point. The function Ghigh has a couple of Yang-Lee branch points on the imaginary axis and is analytic 
at the origin, while the function Giow is expected to have an essential singularity at C = 0. We argue 
that the domain of analyticity of Giow(C) is the wedge | arg((^)| < 7r(l — Ag°^). The scaling functions 
Ghigh and Giow are analytically related at infinity and can be expressed through a third scaling function 
$(r/), where r/ = hT^I'^^ t is the dimensionless temperature coupling. 



^ Even in the Ising case, p = 3, it is not clear how to characterize the unstable phase in terms of the gas of dense 
self-avoiding loops and lines. 

* We find more convenient to introduce two different, although related by a phase factor, dimensionless variables, ^ for 
the HT regime and C for the LT regime. 
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In Sect. 3 we formulate the the analytic properties of the free energy of 0{n) quantum gravity, or 
shortly 0„QG, which describes the continuum limit of the 0(n) model on a dynamical lattice. The 
geometrical expansion in this case involves non-intersecting loops and lines on planar graphs as the 
example in Fig. 6. 0„QG is formally defined as a perturbation of Liouville gravity with matter central 
charge (1.8) by the action 

which generalizes (1.4). The vacuum energy of On QG must be of the form 

' i^^high(e), ^ = t-''h if t>0, 



F{t, h) 



(1.17) 



U-ireiow(C), Q = {-t)-^h if t<0, 
where u is the ratio of the Liouville dressing exponents associated with the two perturbing operators, 

u = — = — - — . (1.18) 
ttg 4 

Again, the functions C/high and C/jow are related analytically at infinity and can be expressed in terms of 
a third scaling function $(??), where r/ = h~^^^t is the dimensionless temperature. This determines 
the form of the expansion at infinity 

oo , oo 

j=0 j=0 

At the origin the function Qhi^ is analytic, while ^low is expected, according to [22], to have a power- 
like singularity. The leading term in the expansion is / ~ ^^/"^s""^ where 0;'°^ is the Liouville exponent 
associated with the magnetic operator in the low-temperature phase. The Z-Z conjecture concerns the 
subleading term, which should be ~ /^'™ , with gio^ given by (1.7). We develop further the arguments 
of [22] by considering the contributions of droplets-within-droplets configurations. In this way we 
were able to predict not only the first two terms, but the form of the whole series expansion the 
scaling function ^low, 

f Slow 

^low iO=Q0 + C,f + C2f (1 + C3 ^ + . . . ) 

_ 4(p-l) 4n (1-20) 

n>0 

In Sect. 4 we compare the speculations of sections 2 and 3 about the behavior of the 0(n) free 
energy with the exact solution, obtained for the microscopic realization of 0„QG as the 0{n) loop 
model on a dynamical lattice. The derivation, based on the correspondence with the 0{n) matrix 
model, is given in Appendix A. The specific free energy J^{t, h) was found in a parametric form, 

yP yP-1 y-(p+l)/2 



^ ' 2 2 1 - i) 1 + i 



- ■ p (1.21) 
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For p = 3 the equation for the 0(n) free energy coincides, up to a shift by a term ~ t^, with Boulatov- 
Kazakov solution of IQG, eq. (1.6). 

From the solution (1.21) we find the scaling functions Qhi^ and ^low They represent two differ- 
ent branches of the same meromorphic function Q, which is defined as the analytic continuation of 
^high(0 under the YL cuts. The function Q{^) has the symmetries 

g{0 = G{-0 = W), (1-22) 

which are inherited by those of the function Qhigh and which preserve by construction the HT sheet. 
Starting from the HT sheet and taking different paths, one can achieve (for general p) four different 
copies of the LT sheet, related by the symmetries (1.22). The HT and the LT sheets are connected by a 
finite number of auxiliary sheets. The extended analyticity assumption, which is fulfilled here, means 
that there is no additional singularities on the coimecting sheets. 

We computed the series expansions of ^low and ^high at C = and at C = oo- The expansion of 
Qiow at the origin is indeed of the form (1.20), which confirms the Z-Z conjecture in its stronger form 
used in Sect. 3 as well as the existence of Langer singularity at the origin associated with the presence 
of a metastable state. If the gravitational field is 'switched off', the power-like singularity at C = 
should turn into an essential singularity associated with this metastable state. Let us stress that we 
expect Langer-type singularity only in the 0(n) loop model which has a representation in terms of a 
gas of non-intersecting polymers. The conventional 0(n) model should not have such a singularity. 

To summarize, in this work we argue that Langer' s singularity in presence of small magnetic field, 
originally observed in the Ising model, is in fact a general feature of the 0(n) loop models. We are 
quite confident in our conclusions when < n < 2, while in the interval — 2 < n < 0, where the 
loop gas does not have statistical interpretation, the situation is less clear. Our exact expression for 
the free energy of the gravitational 0{n) model presents a strong evidence about the existence of a 
metastable state in the low- temperature phase and confirms the Z-Z conjecture about the decay rate 
of the metastable vacuum in presence of gravity [22] in the case when the stable phase has non-zero 
central charge. We found that the scaling function Q{^) for the free energy of the 0{n) loop model on 
a dynamical lattice has a pair of Yang-Lee singularities on the high-temperature sheet, a Langer-type 
singularity of the expected form on the low-temperature sheet, and no additional singularities on the 
connecting sheets. In this sense the scaling function obeys the 'extended analyticity' assumption of 
Fonseca and Zamolodchikov [8]. However, it is not clear if this property will remain true for the 
theory on a flat lattice. This is a very important question which deserves further study. 

2 The vacuum energy of the 0(n) field theory 

The analytic properties of the free energy of of OnFT, except Langer singularity, follow immedi- 
ately from the known scahng behavior near the critical points, where the theory becomes conformal 
invariant, and the assumption that on the first sheet there are no other singularities than the critical 
points. These predictions are in accordance with the numerical results for the cases n = 1 [8, 11] 
and n = [30]. In contrast, as we already discussed, we do not have direct arguments in favor of 
the presence of Langer singularity. This is a conjecture which is justified by the exact solution of the 
model coupled to gravity. 

2.1 Continuous transitions and effective field theory on a regular lattice 

The local fluctuating variable in the 0{n) loop model [24, 25, 31] is an n-component vector S{r) 
with unit norm, associated with the vertices r of the honeycomb lattice, and interacting in an isotropic 
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way along the links <r, r'>. In the presence of a constant magnetic field H, the energy of a spin 
configuration is defined by the "geometric" hamiltonian 

n{S}\ = - J2 log + S{r) ■ Sir')) - log (l + S{r) ■ h) . (2.1) 

<r,r'> r 

The logarithmic form of the nearest neighbor interaction leads to a simple graphical expansion 
which generalizes that of the Ising model and can be mapped to a solvable vertex model. The partition 
function is defined as the trace (the integral over all spin configurations) 

Z^(„,(r;r) = TVe-«[{^>], (2.2) 

where by the 0{n) symmetry Tr:{SaSi,) = 6afi and Ti (S'a) = Tr{SaSbSc) = 0. Expanding the inte- 
grand as a sum of monomials and using the rules above one can write the high temperature expansion 
of the partition function as the grand canonical ensemble of non-intersecting polymers of variable 
length. There are two kinds of polymers: loops with activity n and open lines with activity H^, where 
H = \H\. Apart of the activity, the Boltzmann weight of each polymer is equal to T"^, where L is 
the length, defined as the number of links covered by the polymer. The high temperature expansion 
of the partition function is a triple series in n, 1/T and H^, 

Z^^^^{T,H)= J2 (l/r)^*°*n#i°°P^ii"2#°P^"""«^ (2.3) 

polymers 

The temperature coupling T controls the total length Ltot of the polymers (the number of Unks covered 
by loops or open lines). An example of a polymer configuration is given in Fig. 5. The geometri- 
cal expansion (2.3) of the 0{n) model allows to consider the number of flavors n as a continuous 
parameter. 

When H = 0, the 0{n) loop model is known to have a continuous transition in the window — 2 < 
n < 2, where n can be parametrized by eq. (1 .9). The phase diagram of the loop gas on the honeycomb 
lattice was first estabhshed by Nienhuis [25]. At the critical temperature Tcrit = 2 cos ^^^^^ the loop 
gas model is solvable and described by a CFT with central charge (1.8). This point is usually referred 
as the 'dilute' phase. For T > Tent, the theory has a mass gap. The low-temperature, or 'dense' , phase 
T < Tcrit represents a massless flow toward an attractive fixed point at Tiow = 2 sin . At this 

point the theory is again solvable and described by a CFT with smaller central charge Qqw defined by 
(1.10). 

The spectrum of the degenerate primary fields in the critical and in the low-temperature CFTs is 
given respectively by 

_ (r(p + l)-5p)^-l .iow_ {rp-s{p-l)f-l 

where r and s are integers. The thermal operator e, which counts the total length of the polymers in 

the expansion (2.3), can be identified with the degenerate field $1,3. Added to the action, it generates 
a mass of the loops. The spin operator s coupled to the magnetic field corresponds to a degenerate 
primary field only if p is an odd integer, 

p = 2m-l. (2.5) 

However, it is convenient to extend the Kac parametrization (2.4) also for non-integer r and s. Then 
the magnetic field can be identified with the field of conformal dimensions Ag = Ag = A„ with 
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m = ^{p+ 1). The magnetic operator is the first of an infinite series of L-leg operators ^Lm,Lm {L = 
1,2,...) [32]. Explicitly, the conformal weights of the thermal and spin operators are given by 

A A A A ^ (p + 3)(p^l) 

A£ = Ai3 = — — -, As = Ap+i p+i = — — — . (2.6) 

P + I 2 ' 2 lb P{P + 1) 

Both operators are relevant. The 0{n) field theory, which describes the continuum limit of the Hamil- 
tonian (2.1), is formally defined by the action (1.11), where 

t-T-Tcrit, h^\H\. (2.7) 

parametrize the vicinity of the critical point of the lattice model. When n = I, or p = 3, this action 
coincides with the action (1.1) for the Ising field theory. Depending on the sign of the temperature, 
the perturbation by the thermal operator e generates a flow toward the massive or to the dense phase. 



As we already mentioned in Sect. 1, the low-temperature phase of the 0(n) loop model defined 
by the Hamiltonian (2.1) does not contain Goldstone modes. This is in contrast with the standard 
0{n) model [33, 34], which is defined by the Hamiltonian 

ms}] = E ^(^) • ^(^') - E ^(^) • ^- (2.8) 

<r,r'> r 

In the case of the Ising model (n = 1), the Hamiltonians (2.1) and (2.8) are equivalent, up to 
a redefinition of the two couplings. For general n G [—2, 2], the two Hamiltonians are supposed to 
share the same universal properties in the high temperature phase and at the transition point, while the 
low temperature phases of the two models are believed to be different [28]. For T < Tc, the 0{n) 
symmetry of the model with Hamiltonian (2.8) is spontaneously broken to 0{n — 1), which leads to 
a massless Goldstone phase with central charge Qqw = n — 1. (Since the model is non-unitary, the 
Mermin- Wagner theorem does not apply here.) In contrast, the low-temperature phase of the model 
with Hamiltonian (2.1) exhibits unbroken 0{n) symmetry. 

In terms of the geometrical expansion, the Goldstone phase appears when the loops are allowed 
to intersect. The Hamiltonian (2.1) of the 0(n) loop model is designed in such a way that the inter- 
sections are strictly forbidden. This is why the Goldstone modes do not appear. It is argued in [28] 
that an arbitrarily small perturbation that allows intersections would cause a crossover to the generic 
Goldstone phase. 



2.2 High-temperature regime 

In the high-temperature phase of the 0(n) model, the deformation (1.11) with f >0,h = Odescribes, 
from short to long distance scales, the flow to a massive theory. For finite t, the specific free energy, 
or the vacuum energy density of the 0{n) field theory,^ scales as 

F{t, 0) ~ t^-^' =t^ {h = (},t> 0). (2.9) 

For finite h and t, the specific free energy must have the scaling form 

j>+i 

F{t,h) = t^Gi^^{0, (2.10) 

^If the theory is defined on a cylinder of radius R, the specific free energy per site is determined by the leading asymp- 
totics of the energy of ground state Eo{R), when the radius tends to infinity, F = limij_>oo Eq (R) /R. 
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where the scaling function Ghigh depends on the dimensionless variable 



,._ _ l-Ai _ (3p+l)(5p + 3) 



{ = r"«ft, "0 = ;^= - — - . (2.11) 
When p is an odd integer, as in the case of the DFT, one should take the derivative of (2.10) in p, 

F{t, h) ~ t"' log t + i™Ghigh(0> p = 2m - 1 = odd integer. (2.12) 

Since for a finite positive t the theory has a mass gap, the scaling function should be analytic at 
^ = 0. The high temperature expansion (2.3) involves only even powers of the magnetic field, which 
means that Ghigh is an even function of ^ in the vicinity of ^ = 0, 

Ghigh(e) = Go + + Gii^ + ... for e ^ 0. (2.13) 

On the other hand, the large h finite and t = the theory is again massive with correlation length 
determined by the conformal weight Ag. Hence the asymptotic behavior of the scaling function at 
infinity is 

Ghigh(e) - + ■■■ for C ^ oo. (2.14) 

The function Ghigh(0 can be analytically continued for complex values of ^. Obviously the scal- 
ing function has branch points at some finite ^, otherwise the two asymptotics (2.13) and (2.14) would 
be incompatible. As argued in [34], the general 0{n) model exhibits in the HR regime the same 
analytical properties than the Ising model: the zeros of the partition are situated along the imaginary 
axis, and an infinite strip centered around the real axis is free of zeros. Then, the free energy would 
have two symmetric branch cuts on the imaginary axis, extending from iti^c to itioo, as is shown in 
Fig. 1. The branch points correspond to the two Yang-Lee edges, and in the vicinity of ^ = zti^c 
the 0{n) field theory is expected to be described by the minimal CFT M.2,b with a central charge 
c^^ = —22/5 [4]. The Yang-Lee CFT has only one relevant operator $1^2 with the conformal dimen- 
sion /^^2 = —1/5. The dimension of the corresponding coupling constant is 1 — = 6/5 and one 
expects, as in [8], that near the YL branch points Ghigh(0 behaves as 

Ghigh(0 = 9a{0 + 9b{0 [il + ef''' + ---, (2.15) 
with QA and qb being regular functions of ^. 

2.3 Low-temperature regime 

In the LT regime {t < 0) the theory flows toward a CFT with the central charge Qqw given in (1.10) 
[29, 35]. The perturbation of the low-temperature CFT at t ^ —00, /i = is driven by the operators 
e and s with conformal weights 



Alow Alow Alow Alow 1 (P 3)(p+l) 

p—l -^'^ 16 p{p — 1) 



and coupled respectively to — 1/t and h. The energy becomes an irrelevant operator, while the mag- 
netic field remains relevant. A finite perturbation at i — >■ —00 with this operator leads to a massive 
theory with specific free energy 

F{t,h)r~.h'/(^-^'r) (t^-oo). (2.17) 
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For the Ising model {p = 3) the low-temperature phase contains no degrees of freedom, Qow = and 
^low _ Q jjjg fj.gg energy is F ^ h and the magnetization 

s = -dF/dh (2.18) 

is a constant. The dimension of the magnetic operator becomes negative for p < 3. Accordingly, the 
magnetization s vanishes at ^ ^ for p > 3 and diverges as for p < 3. In particular, 

for p = 2, which corresponds to the case of the dense polymers, the specific free energy diverges as 
F h?"^/^^. This scaling behavior has been tested and confirmed numerically by Saleur [30]. In the 
opposite limit, /i oo, the free energy is again given by (2.14). On the other hand, for finite negative 
t and /i = the free energy behaves as 

F{t,h) - (-l/t)V(i-A'r) ^ (_t)^ {h = Q,t< 0). (2.19) 
For finite values of h and t the vacuum energy is of the form 

F(i,/i) = (-t)'^Giow(C), (2.20) 
where Glow is the scaling function in the low-temperature phase and C, is defined as 

C = (-t)-'^o/i. (2.21) 
At small C the scaling function Giow has asymptotic behavior 

Giow(C) ~ Go + Gi C'/^i-^^"") + . . . forC^O, (2.22) 

while at infinity G\ 

ow ~ (^^/(^ As) gy construction the function Giow(C) is real for > 0. The 
constant term in (2.22) is in general different than the one in (2.13). Since the susceptibility and the 
magnetization are positive, we have Go > and Gi < 0. The relation between Glow and Ghigh will 
be obtained in Sect. 2.5. 

2.4 Langer singularity 

In the Ising model (p = 3) the scaling function Giow is expanded in the integer powers of C but 
the power series is asymptotic. The analytical continuation of the free energy from positive C has an 
exponentially small purely imaginary discontinuity on the negative axis, which is explained by the 
presence of a metastable state for small negative C,. In the language of QFT, one could say that the 
theory contains a global resonance state with complex energy, the 'false vacuum'. The imaginary part 
of the specific free energy can then be interpreted as the decay rate per unit time and volume of this 
resonance state [36] [9] [10]. Langer [5,6] developed a systematic approach, based on the earlier 
droplet models, to evaluate the discontinuity 5F = 2i3F of the free energy on the negative axis. The 
latter decays exponentially for small magnetic field, 5F ~ exp(— const/|C|). 

By construction, the free energy and the correlations in the 0{n) loop model are analytic in p. If 
the perturbative expansion of the free energy at p = 3 is asymptotic, it will be asymptotic also in the 
vicinity of p = 3. Therefore it is very likely that a metastable state and Langer type singularity should 
exist also for p ^ 3. This conjecture will be later confirmed by the exact solution on dynamical lattice. 
Unfortunately it is not clear how to characterize geometrically the true and the false vacuum within 
the loop gas representation, even in the Ising case p = 3. 
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In the Ising model the energy gap between the stable and the metastable vacuum is evaluated by 
analytically continuing the free energy to negative values of h. For general p the leading behavior of 
the free energy for small h is no more linear in (, 

F{t, h) = i-tf-^ (Go + GiC''^ + ...), 2u 1/(1 - A'r), (2.23) 

and the metastable state should be reached by analytical continuation to the rays arg(^) = ib7r/2i/ in 
the complex /i-plane. Here the free energy is again real and the gap 

fit, = F{t, e±-/2-/.) - F{t, h) (2.24) 

is positive. Langer's theory [5] predicts a nnon-perturbative, exponentially small in /, imaginary part 
of F, equal to the decay rate T (per unit volume and unit time) of the metastable vacuum, 

F{t, e'^'/'^'^h) - F{t, ~ iV. (2.25) 

To evaluate the leading exponential factor, assume that initially the two-dimensional plane is filled 
with the metastable vacuum with specific free energy F+f and imagine a hypothetical stochastic pro- 
cess which allows the statistical system to develop in time. The decay of the metastable vacuum will 
then occur through spontaneous formation of droplets of the stable phase in the "sea" of metastable 
phase. The Boltzmann weight of a droplet of area A and perimeter L is 

Wa,l ~ e^^-^^, (2.26) 

where s is the surface tension at the boundary of the droplet. The surface tension scales as cr = 

p+i 

(^o{—t)~, where so is a numerical coefficient. Assuming that the drops are approximatively round, 
A = L'^/Att, the integral in L is saturated by the saddle point Lg.p. = 27rcr//. The saddle point 
value gives the size of the 'critical droplet', which can trigger a decay of the metastable state. For 
droplets larger than the critical droplet the bulk energy starts to decrease faster than the surface energy 
grows and the droplet energetically favorable to grow. The leading factor in the decay rate T is the 
Boltzmann weight of the critical droplet. A more refined calculation [37], taking into account the 
fluctuations of the form of the droplets, gives also the pre-exponential factor, 

r ~ / 6-"""^/^. (2.27) 

Therefore, with the assumptions we made, the scaling function Giow(C) for the LT regime develops 
an essential singularity ~ Izp'^e"!^! " at the origin. The function Giow(C) is analytic in the C-plane 
with a cut starting at C = 0, which can be placed along one of the rays arg(^) = ib7r/2z/, where the 
discontinuity of the leading perturbative term vanishes, or along the negative axis, which will be our 
choice in the following. 

2.5 The scaling function $(77) 

The functions Ghigh(0 G\ow{0, although defined for two different phases, must be analytically 
connected in the regime where the magnetic field is sufficiently large compared to the temperature. In 
this regime it is useful to introduce the dimensionless temperature 

rj = ^-V'^o (2.28) 
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27l/a)n 



SHD 
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SHD 
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Figure 7: Principal sheet of the 77-plane for 
the scaling function ^{rj). The LT wedge 
|arg(— 7/)| < 2tt/ujo is the image of the prin- 
cipal sheet of the scaling function Giow(C)- Ths 
HT wedge | aig{ri)\ < -k jlijjQ is the image of the 
right half plane of the principal sheet of Ghigh ) • 
The two sectors SHD separating the LT and the 
HT wedges are the 'shadow domain' where the 
analytic properties of $(77) are not known. In [8] 
it was conjectured (for p = 3) that the function 
$(77) is analytic also in the shadow domain. 



with Wo was defined in (1.15), and a new scaling function $(r/) by 

F{t,K) = $(7?). 



(2.29) 



This new variable is suited to study the analytical property of the free energy in terms of the tempera- 
ture. In particular, it allows the description of the neighborhood of the critical point where t ~ 0. In 
the regime t > and /i > 0, the variable rj is related to ^ by 



and the new scaling function is expressed in terms of Ghigh(C) as 

$(,^) =,7T^Ghigh(l/r?"o), (r?>0). 
Similarly, in the regime t < and h > 

rj = - (-1/^0 

and 

Hv) = i-v)^GioAl/i-rjr^), (7?<0). 



(2.30) 



(2.31) 



(2.32) 



(2.33) 



(When p is an odd integer, the above relations are slightly modified. They are obtained by taking the 
derivative in p and using the definition (2.12).) 

Since for /i 7^ the theory has a mass gap, the free energy F{t, h) should be analytic in the tem- 
perature t in some finite strip around the real axis. Therefore the scaling function $(7?) is represented 
by a power series 

$(r/) = $0 + ^1 + ^2rf + ■■■ (2.34) 

with a finite radius of convergence near 77 ~ 0. According to (2.31) and (2.33), the coefficients <I>fc 
of the series (2.34) determine the large expansion of the scaling functions Ghigh(0 and GiowCC) at 
infinity, 

l-j(l-Ae) 



Giow(C)= C^^(-C 



'1/(^0 ^ 



l-j(l-AE) 



(2.35) 



14 



In the vicinity of ^ = oo, the function Giow(C) can be obtained from Ghigh(0 analytical continua- 
tion C ^ C = e^*'''^"^. 

The domains of analyticity of Ghigh and are mapped to different (and non-overlapping for 
p>2) domains on the principal sheet of the scaling function $(77). The map ^ — > r/ sends the positive 
real axis of the ^-plane to the positive real axis of the r/-plane. The right ^-half-plane, where Ghigh(0 
is analytic, is mapped to the high temperature wedge (HTW) of the ?7-plane, defined as 

HTW: - < argr? < — . (2.36) 
In particular, the LY branch points are mapped to 

r?± = e±^'^/2-o r/e, r?c = ^y""- (2.37) 

One can choose the branch cuts starting at the YL edges to go to infinity along the rays arg r/ = 
ib7r/2a;o- (Note that the cuts in the ry-plane are not the images of the cuts of the ^-plane.) 

On the other hand, the map C ^ V sends the positive real axis in the ^-plane to the negative real 
axis of the ry-plane. The principal sheet of G^°^ in (^-plane is mapped to the low temperature wedge 
(LTW), 

LTW: - — < arg(-r/) < — . (2.38) 

In the ?7-plane the Langer branch cut is resolved and the Langer singularity is sent to infinity. The two 
edges of the cut are mapped to the two rays arg(— ry) = ±7r/c<;o- The principal sheet of the function 
$(77) is shown in Fig. 7. 

As in the Ising model, the function ^{rj) is expected to be analytic in the sectors HTW and LTW 
of the ?7-plane. The two sectors separating HTW and LTW were called in [8] 'shadow domain'. If the 
extended analyticity, claimed in [8], holds for all p>2, which means that $(r/) is analytic also in the 
shadow domain, then the scaling function can be reconstructed unambiguously from its discontinuities 
along the two Yang-Lee cuts on the principal sheet of the 77-plane. 



3 The vacuum energy of the gravitational 0(n) model 

In this section we obtain the expected analytical properties of the free energy of the 0(n) model 
coupled to 2d gravity. We adjust the arguments of the previous section to the situation when gravity 
is switched on, making substantial use of the Z-Z conjecture about the effect of the critical droplet. 



3.1 Continuous transitions and effective field theory on a dynamical lattice 

The 0{n) model on a dynamical lattice has the same continuous transitions as the theory on a flat 
lattice. In the continuum limit the sum over the planar graphs with a given topology is replaced by a 
functional integral with respect to the Riemann metric gabix, x) on a variety with the same topology. 
Here, we will restrict ourselves to the case of the sphere and the disc, where the metric can be put 
in the form gab{x,x) = dabe^^^^'^-' by a coordinate transformation. On the sphere, the integration 
measure with respect to the scale factor cj), known as the Liouville field, is defined by the action [13] 

Aiouv = ^ f {dcjyfSx + n [ e^'I'd^x , (3.1) 
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and the asymptotics at infinity 



(/>(x,x) ~ -Qlog(xx), Q = g + 1. (3.2) 

The Liouville coupUng constant g is determined by the critical fluctuations of the matter field and the 
'cosmological constant' /x controls the global area of the sphere 

A= [e^'^(fx. (3.3) 



When the topology is that of a disk, the integral (3. 1) is restricted to the upper half plane, with an extra 
boundary term /Xg J e'^dx. The boundary coupling constant controls the length of the boundary of 
the disk, 

L = J.U.. (3.4, 
The Liouville action (3.1) leads to a CFT with central charge 

CLiouy = 1 + 6(5+ 1)V5- (3.5) 

Besides the Liouville field there are also Fadeev-Popov ghosts with total central charge Cghosts = —26 
(see the standard reviews [38] and [39]). Due to the general covariance, the total conformal anomaly 

in the theory of Liouville gravity must vanish, c uouv + Cmatter + Cghosts = 0. When the matter field is 
the 0(n) CFT with the central charge (1.8), this constraint determines the Liouville coupling constant 
g in (3.1) as 

g = ip+l)/p. (3.6) 
The observables in Liouville gravity are the integrated local densities 

- J^A e^'^^cPx , (3.7) 

where $a represents a scalar (A = A) matter field, and the vertex operator e^°"^ is the Liouville 
dressing factor, which takes into account the fluctuations of the metric. The Liouville dressing factor 
completes the conformal dimensions of the matter field to (1,1). Only such, marginal, operators are 
allowed in the theory coupled to gravity. The balance of the conformal dimensions gives a quadratic 
relation between the Liouville dressing charge a and the conformal weight of the matter field, known 
as the KPZ relation [14-16] 

a{g + l- a)/Ag = 1 - A . (3.8) 

The 'physical' solution of the quadratic equation (3.8) is the smaller one. 

In the Kac parametrization (2.4) of theconformal weight A = A^^^ of the matter field, the physical 
solution of the quadratic equation (3.8) is given by 



2p 



In particular, the operators 
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represent the Liouville-dressed thermal and magnetic operators with dimensions (2.6), and the corre- 
sponding Liouville exponents are 

1 - 3p+l 
CHe = Q:i,3 = Q!s = Oip+1 P+1 = — • (j-11) 

P 2 ' 2 Ap 

The 0{n) gravity is formally defined by the action (omitting the ghost term) 

AOnQG = Arit + Aiouv -tO,-hOs. (3.12) 

Since all operators in Liouville gravity are marginal, the dimensions of the matter fields cannot be 
extracted from the coordinate dependence of correlation functions of the theory. In fact, the correlation 
functions in Liouville gravity do not depend on the coordinates at all. The dimension of the matter 
component of an operator C> = $ e'^""^' is determined by the response to a rescaling of the metric, that 
is, to a translation of the Liouville field. As a consequence of the KPZ relation (3.8), the exponent a 
is positive if the matter field $ is relevant and negative when the matter field is irrelevant. 

Let Z{t, h, /i) be the partition function of Liouville gravity on the sphere. The role of the RG time 
in Liouville gravity is played by the zero mode of the Liouville field, which is equal to the logarithm of 
the square root of the area (3.3). Therefore sometimes it is instructive to consider instead the partition 
function on the sphere for fixed area, which is given by the inverse Laplace transform with respect to 
the cosmological constant, 

ZA{t,h) = j^^.e^>^ Z{t,h,ii). (3.13) 

Here the contour of integration goes on the right of all singularities of Z{t, h, fx). In the limit of 
infinite area the partition function behaves as 

ZA{t,h) ^ A-y^ e-^^^''^\ (3.14) 

The exponent is by definition the specific free energy for the 0{n) quantum gravity. 

3.2 Analytic properties of the specific free energy in tlie high-temperature regime 

In the previous section, the scaling behavior of the specific free energy of the OnFT near the critical 
points was determined by the dimensions of the two coupling constants, t and h. After coupling to 
gravity, the dimensions of the coupUng constants reflect the response of the corresponding operators 
to translations of the zero mode of the Liouville field. With the convention that the cosmological 
constant jjl has dimension 1, the dimension of the coupling constant for the operator (3.7) is a. 

To establish the scaling properties of the free energy in 0„QG, one can essentially repeats the 
arguments of the previous section, replacing the dimensions 1 — A of the coupling constants by 
the Liouville exponents a. In the high-temperature regime, t > 0, the specific free energy of the 
gravitating 0(n) field should be of the form 

j-(t,^) = W°^ghigh(0, (3.15) 

and the scahng function ^high(C) depends on the dimensionless strength of the magnetic field, defined 
as 

c 3p-|-l 

^ = t h, io = — = — - — . (3.16) 
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For finite positive t, the theory has a mass gap and the scaling function must be even and analytic near 

ahigh(0 = ^0 + 02^ + 04^^ + ..., m < ^c. (3.17) 

The free energy is expected to have two symmetric Yang-Lee branch cuts on the imaginary axis^, 
extending from = ±i^c to ±ioo, with a value different from the flat lattice model one. Near 
the YL singularity, the theory is that of a Liouville gravity with Cmatter = —22/5, or g = 5/2. The 
magnetic field is coupled to the only primary field ^1,2, which have Liouville dressing exponent 
Q;i,2 = 3/2. In such non-unitary theory, the operator that determines the behavior at large distances 
(here large cf)) is the most relevant operator Oi,2, while the perturbation with ■ Oi^i, 6^ = ^ — 
leads to a 'correlation area' Ac ~ ((5^)-°i,2/q:i,i = ((5^)~^/^. Therefore the specific free energy must 
behave for ^ « as 

gMgh(0 ~ reg. + {e + ecf^ + ■■■ (3.18) 
where "reg" denotes regular terms in 6^. 

3.3 The specific free energy in the low-temperature regime and the Langer singularity 

As in the case of rigid geometry, when i — > — 00, the theory flows to the low-temperature CFT in 
which the coupUng constant of the Liouville field is 

9iow=p/{p-l). (3.19) 
The thermal operator becomes irrelevant, but the magnetic operator remains relevant, 

,T,low — ^low _ low _ low _ 3p — 1 

p—l —'— 4(p— 1) 

A finite perturbation at t ^ —00 with the magnetic operator leads to a massive theory with specific 
free energy 

1 / low 4(p-l) 

J^{h) ~ /i^/"- = h^i^, t -00. (3.21) 
For finite values of h and i < the free energy has the form 

J^(i,/i) = (-i)f giow(C)- (3.22) 

where ^low is the scahng function in the low-temperature phase (analytically continued from a real 
positive magnetic field) and ( is the dimensionless variable 

C = (-t)-'-/i. (3.23) 
According to (3.21), the leading behavior at small ( of the scaling function is 

e?iow(C) = Qo + Qi C^"^^ + subleading terms, C ^ 0. (3.24) 
Since the magnetization s = dF /dh is positive, the coefficient Qi must be negative. 



''The YL singularity in the IQG was first discovered in [18] 
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Assuming that there is a metastable state, its energy is obtained by analytic continuation from the 
positive axis of the ^-plane to the ray arg ( = 7ra'°*, where the leading free energy is again real. The 
energy gap between the metastable and the true vacuum is now 

fit, = T{t, e--^/i) - T{t, h) = -2^1 {-tf C'/"'^ > 0. (3.25) 

Let us sketch the semiclassical argument presented by A. and Al. Zamolodchikov in [22]. Assume 
that in the beginning we have an infinite, globally flat space filled with the metastable phase. The 
corresponding classical solution for the metric is </) = 0. Assume also that there is a stable phase 
with negative, compared to that of the metastable phase, energy density — / and that the decay of the 
metastable phase goes through formation of droplets of the stable phase. 

Let us calculate the energy of a configuration with one circular droplet of the stable phase which 
covers the circle |x| < r. Then the total energy of the droplet is 

^droplet = -/ / e2<^(") +a [ e-^^^^dx, (3.26) 

J\x\<r J\x\=r 

where a is the surface tension at the boundary of the droplet. The metric in the presence of the droplet 
is given by the extremum of the Liouville action (3.1) with g and /x replaced by 

9 5iow =P/(P- 1)' 

(3.27) 



— / if |x| < r, 
if \x\ > r. 



The corresponding classical solution is obtained by sewing the solution of the Liouville equation 

ftow d^d^cl^ix, x) = -f e^-^^^'*) (3.28) 

inside the circle with the flat solution cf) = outside. Since the cosmological constant is negative 
inside the droplet, the general solution of (3.28) with radial symmetry describes a metric with constant 
positive curvature: 

e2<^(^) = (3 29) 

This is the metric of a sphere with radius R, which is related to / by 

^ ^ (-3 30) 

This solution should be sewed with the solution cf) = for \x\ > r. This gives the condition = 
2aR, which has two solutions for a, 

R±VW^ 

a± = ^ , (3.31) 

with an obvious geometrical meaning. Imagine that the sphere of radius R is cut into two disks 
along a circle with radius r. Then the two disks correspond to the two possible solutions for the 
parameter a. The smaller disk (a = a_) resembles a cap and the larger disk (a = a+) resembles 
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a bubble. The saddle point for the integral in r and R gives for the energy of the critical droplet 
= fi'iowlog(l + TTcr^ /^(low/)- For gxowf ^ the unstable direction at the double saddle point 
is mostly along the "perimeter growth" direction r, just as in the case of the flat geometry. In the 
opposite limit yiow/ cr^, the growth of the droplet is mostly in the "inflation" direction R. The 
decay rate 



1 + 



■KG 



5low/ 



"Slow 



(3.32) 



matches in the limit yiow/ ^ (strong metastability) the standard droplet model result (2.27), while 



in the limit 5iow/ *C s (weak metastability) it gives a power law, V 



(510 



2 \ Slow 



To summarize, in the limit of weak metastability the critical droplet looks like a bubble connected 
to the rest of the surface by a thin neck, as shown in Fig. 4. In this case the decay rate behaves as J^'™ 
where ^low is defined by (3.19). This phenomenon was called in [22] 'critical swelling'. 

Although the above argument is valid only in the quasi-classical limit Qqw — > —oo, where ^low = 
— ciow/6, it was argued in [22] that the power law persists for finite ciow, and that ^low defined by 
(3.19) gives the exact nucleation exponent. This is supported by the following heuristic argument. In 
the limit of weak metastability the classical solution for the metric inside the critical droplet is that of 
an almost complete sphere. One can speculate that in the quantum case the metric within the critical 
droplet starts to fluctuate, but is still connected to the rest of the surface by a thin 'neck'. If this is 
true, it is easy to evaluate the exact contribution of the droplet to the free energy. In Liouville gravity 
this is the partition function on a sphere with a puncture and a cosmological constant — /, which is 
known [ 14] to be ~ /^°* . This argument is however impeded by the negative sign of the cosmological 
constant. 




Metastable vacuum 




Figure 8: A nested droplets configuration. 



Figure 9: The metric of the nested droplets 
configuration in presence of gravity. 



The Z-Z conjecture formulated above gives a prediction for the subleading term in the small C, 
expansion of ^low, provided the leading term is known. From (3.25) we find 



aiow = ao +Ci/ + C2r™ + ... 

^ ^ 4(p— 1) ^ 4p 



(3.33) 



with ^ < 0. We can develop further the arguments of [22] by taking into account the droplet-inside- 
droplet configurations as the one shown in Fig. 8. The solution of the Liouville equation for such 
a nested configuration describes a metric that resembles the surface of a cactus (without the spines), 
depicted in Fig. 9. It consists of several bubbles of the stable phase connected by thin throats of the 
metastable phase. The contribution of n droplets within a droplet factorizes into the partition function 
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of a sphere with n + 1 punctures and n partition functions on a sphere with one puncture. The total 
power of 1^1 is 

JSlow-n^y9lo„)n ^ ^3^ + 3^^ n = 1, 2, . . . . (3.34) 

In this way the Z-Z conjecture gives not only the subleading term, but also the complete small 
expansion of the scaling function Q\ov/- 



aiow(c) = Oo + J2 ^n+i c'^^^ . (3.35) 



^ 4(p-l) I 4n 

n>0 

vlow _ „ 3p-l 



The position of the Langer cut is at arg^ = ira^ = tt ^^Ziy When n > 1 the Langer singularity 
appears before reaching the negative axis, while for < n < 1 it appears after crossing the negative 
axis. 

3.4 The scaling function $(17) 

To study the free energy for a large magnetic field, or equivalently at a small temperature t, we 
introduce the scaling function $(77) defined by 

J^{t, h) = /i^/"=$(?7), ri = h-^'"^ t . (3.36) 

Depending on the sign of t, the new scaling function $(77) is related to ^high(0 or ^iow(C)' 

mn) = ["^"^^ ^h:gh(l/rr) if i > and > 0, 

if t <Oand/i>0. 

As in the flat case, for a given value /i / of the magnetic field, the matter theory has a mass gap 
and Fit, h) should be analytic in t. Hence, the scaling function <I>(r/) should be analytic in a finite 
strip containing the real axis. Near 77 = 0, it is represented by the power series 

$(r/) = $o + $i?? + $2r?^ + .-- (3.38) 

with finite radius of convergence. According to (3.37), the series expansion of ^{rf) near the origin 
determines the expansion of ^high(0 ^iow(C) ^f infinity: 

00 , 00 

^high(o = Y.^3 ^iow(c) = c"^- (3.39) 
3=0 3=0 

The principal sheet of $(r/) can be split into four sectors as in Fig. 10. The right half -plane 
> 0, where ^high(0 is analytic, is mapped to the the HT wedge on the principal sheet of ^{rj), 
defined as 

HTW: I arg??| < — = -— . (3.40) 
2lo 3^? + 1 

The two Yang-Lee branch points in the ^ plane correspond to two simple branch points of $(77) at 



V± = e^'^'-^'^rjc, Vc = ic'""- (3.41) 
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Figure 10: Principal sheet of the ?7-plane for 
the scaling function ^(r/) in presence of grav- 
ity. The LT wedge |arg(-r/)| < 3^^^ is 



the image of the domain of the analyticity of 
the scaling function C/iow(C)- The HT wedge 
I arg(7])| < is the image of the right half 
plane of the principal sheet of tJhigh(0- 



To exploit the analyticity of ^*(r/) at sufficiently small ry we place the two YL cuts such that they go to 
infinity along the rays arg(?7) = ib7r/2a;. Furthermore, the ^-plane cut along the negative axis, where 
^iow(C) is analytic, maps onto the LT wedge on the principal sheet of <I>(r/) defined as 

vr Ait 

LTW: I arg(-7?)| < - = — — . (3.42) 

u 3p + 1 

The two edges of the Langer cut in the (^-plane are mapped to the two rays arg(— t/) = ztTr/w. 

Coupling to gravity, or taking the average with respect to geometries, cannot add new singularities 
to the free energy. Therefore, if the OnFT satisfies the extended analyticity, this must be true also for 
OnQG. If this is the case, the function must be analytic in the whole 77 plane with two cuts going 
from the points (3.41) to infinity, including the 'shadow domain' which separates the HT and the LT 
wedges. 

4 Analysis of the matrix model explicit solution 

In this section we show that the free energy extracted from the exact solution of the discrete model 
does have the properties claimed in the previous section. 

4.1 Lattice formulation on a planar graph 

The 0{n) loop model can be defined on any trivalent planar graph T. For any value of n, the partition 
function is equal to the sum of all configurations of self- and mutually avoiding loops with activity n 
and open lines with activity H"^, 

Z (T H'V) ^ ^ i^-j^y^^yotal length of the polymers ^#loops open lines 

polymers on F 

An example of such a configuration on a graph with the topology of the disk is given in Fig. 6. 

In the 0(n) model on a dynamical lattice, the planar graph F itself is treated as a statistical object. 
The partition function of the 0(n) model on a dynamical sphere of volume N is defined as the average 
of the partition function (4. 1) over all graphs Fjv having N vertices and the topology of a sphere, 

2-sphere^j.^ ^) = X] ^O(n) ^5 ^ n) ■ (•4 2) 
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Similarly, the partition function on a dynamical disk of volume N and boundary length L is defined 
as a sum over all planar graphs F jv,l with N vertices and L external lines, having the topology of a 
disk, 

^N%{T, H) =^ (T, H; Tjv.l) • (4 3-) 



Tat, 



L 



The disk partition function above is defined for the simplest boundary condition with no special weight 
assigned to the boundary spins. This boundary condition means that the loops are repelled from the 
boundary, as in Fig. 6. 

We are going to compute the universal piece of the specific free energy 

T{T,H)':^' lim ^-^^'ty^H) ^ 

as well as the boundary specific free energy 

- def Z%'l{T,H)-NT(T,H) 
F^{T,H)= lim ^'^^ ' (4.5) 

N,L-^oo L 

For that it is easier first to evaluate the grand partition functions 

00 

N=0 

00 CX) 

JV=0 L=0 



(4.6) 



where /2 and fl^ are the lattice bulk and boundary cosmological constants. The large N and L behavior 
of the micro-canonical partition functions (4.2) and (4.3) is determined by the behavior of the grand 
partition functions near the rightmost singularity in p, and /Zg, which we denote respectively by jl'^ and 

The universal free energies are determined by the behavior of the partition functions in the vicinity 
of the critical point T = T'^, H = 0, which is parametrized in terms of the renormalized couphngs 

fi^fl-fi'', li^^K-Ul t^T-T'^, h^\H\. (4.7) 

All the information about the universal behavior is contained in the singular parts of the partition 
functions, which we denote by Z and U, 

^sphere(-^ T, H) = Z^i^^^^p, T, H) + Z{^i, t, h), 

Z^''\jl, fl^,T, H) = Z^^ifi, jl^,T, H) + U{^l, l^^,t,h). 

The regular terms describe degenerate configurations that diverge in the large N limit and should be 
thrown away in the renormalization process. The renormalized coupling constants t, h, fj,, ji^ can be 
identified with coupling constants of Sect. 3 up to numerical normalization factors. 
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4.2 The solution 

We give the derivation of the analytic expressions for Z and lA in appendix A. The solution for the 
disk grand partition function lA is analytic in the boundary coupling everywhere in the complex 
plane, except for a branch cut along the interval [—00, —M\. The function M = M(;U, K) gives 
the position of the rightmost singularity of U in the [i^ complex plane. The renormalized boundary 
cosmological constant /Xg is defined so that at the critical point |Lt = t = the singularity is placed at 
the origin: M(0, 0, 0) = 0. 

The solution can be written conveniently in a hyperbolic parametrization which resolves the 
branch point at [i^ = —M. The boundary coupling //g is replaced by a hyperbolic parameter r, 
defined as 



/Xb =Mcosh(T). (4.9) 

All the observables must be even functions of r, and their expression can be written solely in terms of 
hyperboUc cosines. To shorten the formulas, we introduce the notation 

Ca(r) = ^cosh(ar), a G M. (4.10) 

Then, up to a normalization of the coupling constants and the partition functions/ the function 
M(/x, t, h) is determined by the transcendental equation 



2 S(i-i)+ i + i 

and the derivatives of the disk partition function in /Xg and in ji are written in the parametric form as 



l^=^-t—. TT + h'^ 1— > (4-11) 

^2 2(1 -i) 1 + i 



r,j, M "p^ Cp-1 (r) - M Cp+i (r) 

-8K ^ ^-(-) + ^^-(-) + ^YpiMP 

dU (4-12) 

Mb = Ci{t). 

Finally, the second derivative of the partition function on the sphere is given by 

u = -d^Z = M^/P. (4.13) 

One can eliminate M from (4.11) and (4.13) and write the equation of state for the susceptibility 
u = -dlZ, 



yP yP-l -(p+l)/2 

u = t — T- + ^ . (4.14) 

V p' p 



' The signs cannot be changed by a renormalization. They follow from the requirement that the divergent singular parts 
{—d^^Z, {—dx)^U and {—8^)^14 are positive. (The singular parts that do not diverge contain a finite regular positive part 
and do not need to be themselves positive.) 
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4.3 Microcanonical partition functions and analytic expressions for the bulk and 
boundary specific free energies 

The specific free energy can be extracted from tlie leading exponential beliavior of tlie micro-canonical 
partition function when the volume tends to infinity. Knowing the grand canonical partition functions 
in the thermodynamical limit, we can reconstruct the micro-canonical ones. In this limit, the discrete 
sums in (4.6) are replaced by integrals with respect to the area A and the length £: 



poo 

Z{ii,t,h)= dAe-^""^ ZAit,h) 

° (4 15) 

U{ii,fi^,t,h) = dA die-"^-!"^^ ZA,e{t,h). 
Jo Jo 

The partition functions for a given area and length are computed by taking the inverse Laplace trans- 
forms. With this definition the regular part of the specific free energy is automatically subtracted. 
For the partition function on the sphere, Za, we find, after integrating by parts. 



ZA{t,h) = j^ l ^^^'^dlZ{t,,,,h) 



where ix{u) is given by (4.14). The integral over ji goes along a contour going upwards and having all 
the singularities of on the left. The large A asymptotics, given by the contribution of the saddle 
point, is 

ZA{t,h)r. ^-3-1/2 g-.F(t,ft)A^ 

where J^{t, h) is by definition the specific free energy, evaluated by the saddle point condition 

J^{t, h) = -^liy), {dp./du)u=y = 0. (4.18) 

The exponent for the power in the saddle-point result (4.17) corresponds to a pure gravity theory 
(cmatter = 0). This is because the asymptotic is evaluated for values of the area and the length much 
larger than the correlation volume of the 0{n) spins, A^^or = h). 

In a similar way, the micro-canonical partition function on the disk Za/ is given by a double 
integral in /i and . The large i asymptotics is determined by the rightmost singularity of U in the 
//g plane, which is the branch point at /Ug = —M. This gives for the boundary free energy 

Mt, h) = {dfi/du)u=y = 0. (4.19) 

The bulk and the boundary specific free energies are intensive characteristics of the system and as 
such do not depend on the global properties of the random surface such as the number of handles or 
boundaries. 

From the equation of state (4.14) one obtains the analytic expressions for the specific free energies, 

2 2(1 -i) l + l 

__l p-i f 3p + l _ pi3p-l) ;\ (4.20) 
" 2^ [p+l^ p^-1 )' 

where the saddle point value u = yis the solution of 

y(Sp-^)/^[y- t)=h\ (4.21) 
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4.4 The scaling functions ^high(£) and ^iow(C) 

Let us first analyze the scaling function for the HT regime (t > 0), defined by (3.15), 

T{t, h) = tP ghigh(0, e = h t-^^P+^^/\ (4.22) 
By (4.20) and (4.21) the scaUng function is given in a parametric form by 

The solution can be expanded (see Appendix B) in a Taylor series in 



(4.23) 



3p+l, 



2^r[3^n + 2-p] ^' 2(p-l) 



The expansion is convergent in the circle |^| < ^c, where 

(3p+l)^ 

is determined by the position of the nearest singularity d^/dy = 0. The function y(^) has in general 
an infinite number of branch points, all at distance from the origin, 

6,± = ±ie^'^'=^ec, keZ. (4.26) 

The two branch points on the physical sheet, = ^o,± = ^i^c, are of course the positions of the 
two YL singularities. Near the Yang-Lee branch points the scaling function has the expected behavior 
of the YL model coupled to gravity, eq. (3.18). 

Looking at the integral representation (4.16), it is easy to see that the YL edges appear as the result 
of condensation of zeroes of the microcanonical partition function ZA{t, h). When ^ is close to 
or to , the second derivative vanishes at the saddle point together with /i'(n) and the integral 

is approximated by Airy integral. When |^| < the fixed area partition function has no zeroes and 
the power 3/2 singularity of ^high(C) reflects the asymptotics of the Airy function for large positive 
values of the argument. On the other hand, when ^ is purely imaginary and |^| > the product of 

YL 
+ 



the Airy functions associated with and has zeroes with density 



P^A^^^fe + Wc, (4.27) 

which condense into the two YL cuts. 

We saw that the scaling function ^high(0 the analytic properties anticipated in Sect. 3.2, 
for which there was little doubt. Now let us turn to the scaling function ^iow(C)' whose analytical 
properties were the subject of the speculations of Sect. 3.3. 

The scaUng function ^iow(C) for the LT regime {t < 0), defined by (3. 15), 

Ht,h) = {-trgioAO, c = h{-t)-^^p+'^/\ (4.28) 
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is written in parametric form as 



1 f p{3p - 1) 3p + l 



(4.29) 



y 



(3p-l)/2 



+ 



The asymptotics for small and large positive ( have the form anticipated in Sect. 3.3, 



^low(C) = < 



4p 



C^p+i if C ^ oo, 

4(p-l) 
C 3p-l if ^ ^ . 



(4.30) 



On the principal sheet, the cut relating the branch points at ^ = and ^ = oo is placed on the negative 
real axis. The small ( expansion of the bulk free energy is computed in Appendix B. It is indeed of 
the form (3.35), 



+ 



P 



p(3p- 1) J(P-i) 1 Mp-1) 4 



4(P-1) 8_ 



(4.31) 



The coefficient Qn of the series are given by 



^0 = 0; Qn 



1 p 



r 



n! 3p — 1 Y 



3p+l _ P+1 
Sp-l"" 3p-l 



(n > 1). 



(4.32) 



The vanishing constant term reflects our convention for the integration constant in the derivation of 
the transcendental equation for /x in Appendix A.2. A non-zero constant term is generated by a 
redefinition of the cosmological constant fi ^ iJ. + c\t\P. 



4.5 The scaling function $(77) 

Here we analyze the scaling function $(77) defined by (3.36), 

J='{t,h) = ^{r]), r] = tK'^. 

Using (4.20) and (4.21) we write the parametric representation of $(r/). 



2 p^-i 

Zp-l 

rj = V — V 2 . 



2 p-i 



The function $(77) is analytic inside the disk 



where it is given by the power series (3.38). For the coefficients we find (Appendix B) 



3p + l 
"2(p + l) 



p i-iy 



r 



2n 2p_ 

3p+l 3p+l 



Sp+1 n! Y 



Ap+2 _ 3p-l 
3p+l 3p+l 



n 



(4.33) 



(4.34) 



(4.35) 



(n > 1). (4.36) 
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It/5 3;i/5 5^/5 7^/5 9it/5 2n 

Figure 11: A circular slice of the Riemann surface of $(77) for to = 2 (p = 3). 



For general p the function ^{t]) has an infinite number of simple branch points, the solutions of the 
equation dv/drj = 0, 



Vj = e^^'^' Vc, j e (4.37) 

There are no other singularities at finite t]. The cuts starting at rjj go along the radial direction and 
end at the branch point at r/ = 00, which is in general of infinite order. The principal sheet has two 
branch points, r]_^_ = r]Q and r]_ =771, which are the images of the YL branch points in the ^ plane, 
as shown in Fig. 10. The two branch cuts separate the sectors HTW and LTWUSHD defined in Sect. 
3.4, in which the scaling function has different asymptotics, 

( -riP if ry ^ 00, 77 G HTW, 

^'(7?) ~ <^ ' p+i (4.38) 
[-(-r/)3^ if 77 ^ 00, 77 G LTWUSHD. 

According to (4.34), the Riemann surface of $(77) is symmetric under a rotation at angle g^^- 
Therefore all sheets are rotated images of the principal sheet. Let us label the sheets by the integers 
. . . , S-i, So, Si, . . . , where 5*0 is the principal sheet. Then the sheet 5^ has two cuts going in the 
radial direction and starting at the points rjk and 77fc_|_i. For irrational p, the Riemann surface has an 
infinite number of auxiliary sheets. 

Let us now consider the case when p = 2m — 1 is odd integer, when the matter QFT is a pertur- 
bation of a unitary CFT The Riemann surface of $(7/) has '^^^^ = 3m — 1 sheets and 3m — 1 simple 
branch points 770, ... , ri3m-2- At infinity, the 3m — 1 cuts join in a ramification point of order 3m — 2 
where ^> ~ ^m/{3m-2) sheets along a circle with radius larger 7/c, the exterior part of the 

Riemann surface will split into 2 disconnected pieces as shown in Fig. 1 1. When analytically continu- 
ing from the HTW, one returns to the principal sheet after circling once the origin. When analytically 
continuing from the LTW, one should circle 3m — 2 times around the origin in order to return to the 
principal sheet. 

4.6 The complex curve of the specific free energy 

The two scaling functions were obtained by eliminating y from the two equations (4.20) and (4.21), 
assuming that the temperature coupling is positive or negative. The equations allow to consider the 
coupling constants t and h as non-restricted complex variables and define the scaling functions in the 
HT and the LT regimes as the restrictions to two different sheets of the same analytic function G{^), 
defined by (4.23), 

(J _ 1 P(3p-1) p-l _ 1 3£+l p 

^ ~ 2 y 2 p+1 y .4 
e = y^'^-'^" (y-i). 



2j-l 
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In this sense one can speak of the complex curve of the specific free energy. When p and 2ijj = 
are integers, this is an algebraic curve of genus zero. 

The scaling function t/high(0 the restriction of G{0 the principal (HT) sheet, while the 
scaling function is obtained as the analytic continuation of ^high(C) below the YL cuts, 

^high(0 =^(0, HT sheet, 

^. 3P+1 ^. 3p+i (4.40) 

QioAC) = e^*"^ ahigh(e^*"— C), C = e^*"— e G LT sheet. 

Let us first remind the structure of the Riemann surface of Q{C) in the simplest case of the Ising 
model, which has been analyzed in the unpublished work by Al. Zamolodchikov [21]. When p = 3, 
the Riemann surface represents a five-sheet cover of the complex ^-plane. The five sheets of the 
Riemann surface of the curve C'^ = — are shown in Fig. 12. There are two copies of the LT 
sheet, each of then connected to the HT sheet by an auxiliary sheet having two cuts. After changing 
the variable ^ ^ C, = e^*^'^/^^, the Langer cut in the two copies of the LT sheet places itself on the 
negative axis in the ^-plane. 



© 
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Figure 12: The five sheets of the Riemann surface for the Ising model, p = 3. 

For general p the analytic continuation from the HT sheet can be performed in four different ways, 
by going in clockwise/anticlockwise direction beneath the upper/lower YL cut. That is, the HT sheet 
is connected in general to four different copies of the LT sheet, which form an orbit of the Z2 x Z2 
symmetry (1.22) of the meromorphic function Q{i). This symmetry preserves the HT sheet (or each 
of the copies of the HT sheet, if there are several of them), but not the LT sheets. 





Figure 13: Analytic continuation under the 
YL cut from the HT sheet to the LT sheet. 



Figure 14: Analytic continuation under the 
YL cut from the HT sheet to the LT sheet. 



In Fig. 13 we show the part of the Riemann surface containing a path, starting from the real axis on 
the HT sheet, going under the upper YL cut and continuing anti-clockwise, until it reaches the positive 

• 3p+l 

axis {C, = e 4 ^ > 0) of the LT sheet. For convenience we place the Langer cut on the negative 
(■-axis of the LT sheet. The total angle swept by the path is a;7r = ^^^^vr. We taylor the sheets of the 
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Riemann surface visited by the path according to the following decomposition of the angular distance 
along the path, 



3p+l 



TT 



Ip + UpTT + TT . 



{up integer, < 0p < vr) 



(4.41) 



This means that the second sheet has, besides the YL cut, a cut along the ray arg(^) = ^vr + 6p, and 
possibly other cuts. The remaining Up identical sheets have two cuts along the rays arg(^) = ^tt + 9p 
and arg(^) = ^tt + 9p + tt. The axis > of the LT sheet is rotated with respect to the positive axis 
of the ^-plane at angle | + 0p if Up is odd, and | + 6*^ + vr if np is even. The image of the path in 
the r^-plane starts at the positive axis and ends at the negative axis, staying all the time in the upper 
half of the principal sheet of the function <&(r/), as shown in Fig. 14. If the path is continued further 
anti-clockwise (the punctured line), it will visit once again the Up sheets and then the second sheet, to 
reach (in general another copy of) the HT sheet. The image of the path in the //-plane returns to the 
positive axis. 

If one continues in the same way, the same pattern will reappear, rotated at angle 29p, and so on. 
The symmetry ^ — )• e^*^p^ of the Riemann surface of Q{^) reflects the symmetry ry — )• e'^^rj of the 
Riemann surface of ^{r]). For rational values of p, the path will eventually return to the starting point 
after a finite number of turns around the origin. 

Let us consider some particular cases. 



• Polymers (n = 0, p = 2). 
The partition function is that of a grand canonical ensemble of open linear polymers. The Riemann 



surface for curve ^'^ = y^/'^ — y^/"^ has 7 sheets, depicted in Fig. 15. There are two copies of the HT 
sheet and four copies of the LT sheet. There is only one connecting sheet on which the function y{^) 
has four simple branch points at it^ci i^^o as well as at ^ = and ^ = oo. 



HT 


















HT 













Figure 1 5 : The seven sheets of the Riemann 
surface for linear polymers, p = 2 (up — 0). 



• p = 2m — 1, m G N. 

In this case the critical theory is a unitary CFT and the operator coupled to the magnetic field is 
the primary field ^m,m- When p is odd integer, the equations (4.39) become algebraic. The function 
Q{^) has, apart of ^ = and ^ = oo, only two branch points, ^ = ±i^c- The Riemann surface has 
3m — 1 sheets, which include one HT sheet and two copies of the LT sheet. The 3m — 1 sheets of 
the Riemann surface are shown in Fig. 16, for m even, and in Fig. 17, for m odd. In the first case 
rip = |m — 3 and 9p = tt, while in the second case Up = and 9p = ^tt. 
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Figure 16: The Riemann surface for p — Figure 17: Riemann surface for p ~ 2m— 1 

2m — 1 with m even. Here np = |m — 3, with m odd. Here Up = 6p — ^tt. 

0p^n. 
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A Matrix model solution for 0(n) spin system in magnetic field 



A.l Mapping to a matrix model and loop equation for the resolvent 



The 0{n) matrix model [23], with 
an extra term taking into account 
the magnetic field, is the zero- 
dimensional planar field theory whose 
Feynman graphs are built according to 
the Feynman rules shown in Fig.(18). 
The partition function of the matrix 
model represents the integral over one 
hermitian matrix M and a 0{n) vector 
Y whose components are hermitian matrices. 




1/J3 \IP P PT ^Ha 

Figure 18: Feynman rules for the 0{n) matrix model in the pres- 
ence of a magnetic field. The magnetic term H ■ Y can be seen as 
a source for tadpoles. 



(A.l) 



At this stage the number n is considered integer. Expanding the free energy of the matrix model as a 
sum over Feynman diagrams, we obtain a sum over all polymer configurations on the disk as the one 
shown in figure 6. The weight of these diagrams is given by 



N 



2~2g 



-#vertices/2 



rp-Ltot j^#loops ^2#open lines 



(A.2) 
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where N is the size of the matrices and g is the genus of the 'fat' planar graph. The functions defined 
in (4.6) are expressed in terms of the matrix model as 

^^P''-^(/z,r,i/)= lim ^logZ^,, 

""'^^^ ^ (A.3) 
Z'^'^fi, ,T,H) = ^Hrn^ _ ( IV log (l - g-^B M) > , /3/iV = e^^ fixed. 

It is useful to shift the matrix variable as M = X + ^I. Integrating over Y, we write the partition 
function (A.l) as an integral over the eigenvalues of X, 

^ f fr -igi- e-/^^("') TT ~ . (A.4) 

with (after rescaling of H) 

F(x) = --(x+l)+5(.+ i) +- (A.5) 
The basic observable in the matrix theory, the resolvent 

is equal, by (A.3), to the derivative of the disk partition function with respect to the boundary cosmo- 
logical constant. One can check that the critical value of the later is e^B = 1/2, so that the spectral 
parameter x in (A.6) is equal, up to a normalization, to the renormalized boundary cosmological 
constant, x = /X3 . 

The resolvent W{x) splits into a singular part w{x) having a branch cut on the real axis, and an 
analytic part Wreg(a;): 

W{x) = Wreg(a;) + W{x) . (A.7) 

The singular part of W is proportional to the derivative of the disk free energy in the continuum hmit 
U{ii^), defined (4.8), 

w{x) = -U'{x). (A.8) 

The regular part is given by 

2V' ix) - nV i-x) 2 ..r.. 

^reg(a^) = J3^^2 = - i2 + n)x-' -ao-aix-a2X . (A.9) 

The function w{x) satisfies a quadratic identity that can be obtained using either the loop equation 
technique or the saddle point approach, 

w^(x) + w'^(—x) + nw(x)w(—x) A 9 9 A 

\ ^ „ ^ = C_4 + C_2 X-^ + Co + C2X^ + C4X^ . (A.IO) 

4 — 
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The coefficients Cfc are functions of T, fl,H, Wi = I3~^ {trX), = I3~^ {tiX-'^) and W-2 = 
/3~^(trX~^) and can be evaluated from the large-x and small-x asymptotics 

2 Wi II- _3, 

wix) = a2X +aix + ao-\ 1 ^ + i- + C(x ), x ^ oo 

X (2 + njx^ 

w(x) = ^ -W-i-W-2X + aix + 0(x^), x^O (A.ll) 
x"^ 



in (A. 10). Obviously 



v2 17-2 

C-4=^, bo=ir—. (A.12) 
2 — n 2 + n 



The solution of the loop equation (A. 10) with the asymptotic (A.l 1) is given by a meromorphic func- 
tion with a single cut [a, h] on the first sheet, with a < 6 < 0. The coefficient c_4 is proportional 
to if^, C2 depends only on T and C4 is a constant. The coefficients, c_2 and cq depend on T, Jl^H 
directly and through the unknown quantities W\, W-i and W-2- Fortunately, we do not need to com- 
pute these coefficients explicitly in order to find the solution in the scaling limit, as we will see in the 
next section. 

A.2 Solution in the scaling limit 

We will solve the loop equation in the vicinity of the double singularity T = T^nt, H = and 
X = 0, Jx = Ji'^. This vicinity is parametrized by the renormalized couplings t ^ T — Tent, h ^ H 
and fiB ^ X, ^ p. — flc- We will use the parametrization (1.9) of the dimensionality n, with 

1/p = 6: 

n = 2cos(7r0), 9 = 1/p. (A.13) 

We will choose a specific normalization of the solution which simplify the formulas. We normalize t 
and h in such a way that the constants 60 and C2 are given by 

Here we used the fact that C2 depends only on T.Once this normalization is fixed, the quadratic identity 
(A. 10) determines the solution uniquely. 

In the continuum Umit, the left bound of the interval supporting the eigenvalue density is sent 
to —00, and we can use an hyperbolic parametrization to resolve the branch cut of the resolvent on 

[-00, -M], 

//g = X = Mcoshr . (A.15) 
The solution will be expressed in terms of the functions 

Cair) = cosh(ar), a = k±9, k = odd integer, (A.16) 

a 

which form a complete set in the space of functions having the required analytic properties. In partic- 
ular they satisfy the shift equation 

a(T + ivr) + a(r - iir) + nC^ir) = (A.17) 
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also obeyed by the resolvent w{t). Note also that the parametrization (A. 15) can be written as x = 
Ci(t). In terms of the parameter r, the functional equation (A.IO) becomes a shift equation for the 
entire function w{t) = w[ijLb{t)]: 

w'^(t + iir) + w^(t) + nw{T + m)w{T) c-^ c_2 , „ , „ ^,^2 , . .c.. 
A • 2, a\ = \n i^w^ + \n (^W2 + co + C2 Ci(r) . (A.18) 

It is noted that the term proportional to C4 is dropped in the renormalization process. 

From the scaling dimension of the resolvent, w ^(1+^)/^ (by definition /i has scaling dimension 
1) it follows that the r.h.s scales as [i^^^ . From here, one finds the scaling of the coefficients: 

hr^H^, tr^n^, co~M^+^ c_2~At^'^^. (A.19) 
We will try to satisfy the functional equation (A.18) by the following ansatz, 

w{t) = Ci+e(r) + tCi_e(r) + /i (A.20) 

The coefficients b± are fixed by the two leading terms in the small x behavior (A.ll), which can be 
transposed into conditions on w{t) near the point r = ±i7r/2, 

6± = ^1m-(^±^). (A.21) 
For this solution the coefficients c_2 and cq are 



h 

C-2 = - 



4(1-^2) 



h^M-^ + M^+^ t M^- 

1 -\- u 1 — a 



(A.22) 



where M = M{fx,t,h). 

We cannot determine the function M{ijL,t,h) by comparing (A.22) with the r.h.s. of (A.IO), 
because the information about the behavior at large x is lost after taking the continuum limit. Instead, 
we will make use of the fact that the derivative dfj_w{x) does not depend on the potential and have, up 
to a normalization, a standard form [40] 

5,^x) = -M-+^^. (A.23) 

smnr 

Comparing this expression with the derivative of the solution (A.20) with respect to jJL at fixed x, 

dM- = -9^cM (m^ - tM-' - M-3) ^lllj^, (A.24) 

yields a compatibihty condition, which can be written, after integration, as a transcendental equation 
forM: 

The cosmological constant is determined up to an arbitrary integration constant ~ t^^^ . 
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The loop mass M is proportional to the second second derivative u = —d'j^Z of the partition 
function on the sphere in the continuum limit, defined in (4.8). There is a simple heuristic argument 
to see that. The derivative —dijJA{iJL^) is the partition function on the disk with one marked point in 
the bulk. In the limit /Xg — > oo the boundary length vanishes and the boundary can be replaced by a 
point. The leading term in the limit /Xg ^ oo is therefore proportional to the second derivative d'^Z. 
Expanding at /xb — > oo we find 

- d^M ~ -jil + M^^ fj,-^ + lower powers of (A.26) 

(the numerical coefficients are omitted). We conclude that the second derivative of the sphere partition 
function is 

u=-dlZ = M^^. (A.27) 

Now we can write the equation of state for the 0{n) model on the sphere in its final form 

1 i-e i+e 

^=^-t^r^,+h'T^, u = -dlZ. (A.28) 



Of course, this equation implies certain normalization of ji and /x 



B" 



B Series expansions of the scaling functions 



We have to deal with the transcendental equations of the form 

H = z + \z'^, a>l. (B.l) 
The solution of this equation as a power series can be derived using the Lagrange inversion theorem. 



^' n\ r[(a-l)n+'7 + l] ^' 

n=0 '■^ / 1 J 



(-irA" r[«n + 7] Ma-i) 



(B.2) 



HT expansion: We would hke to find the series expansion in ^ of ^high, defined by the eq. (4.23). 
It is easier to first consider the derivative 



P 



y-(P+l)/2^ ^2 ^ y(3p+l)/2 _ y(3p-l)/2_ 



de P+l 

After changing the variable to y = 1 + /, we write these equations as 



V P+l 4 2 3p+l 

^ ySp-l^ ^ip-^ = f3p-l _|_j3p-l_ 



p + l 



(B.3) 



(B.4) 



We can now use the result (B.2) with fj, = ^p-i , X = 1, z = /^p-i , a = ^^^1 and 7 = 2±i to write 



(B.5) 



dQtugh p Y 




i-er 


de ~ 2^ 

n= 




n! 



Finally, integrating once with ^high(O) = l/(2(p — 1)) we get 



^high(0 = ^ X] 



r 


p] 


i-er 


r + 


p 





(B.6) 
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LT expansion: For the derivative of the scaUng function ^iow(C)' defined by (4.29), we find 



dGio 



— y 2 . 



P + i 



e=y 



(3p-l)/2 



+ • 



(B.7) 



The equation for ( is equivalent to (B.l), with X = 1, jj = C,"^, z = y ^2 , a = l^rj, 7 = "^zy^ 
by (B.2) we get, after integrating over 



^low(C) 



P 



3p 



n=0 



r 


3p-l" 3p-l 






3p-l" ' 3p-l 



4(p-l) I 4w 
^ 3p-l """Sp-l , 



(B.8) 



The other expansions for the sheets having a branch point at C = are obtained by considering the 
other roots of the function ^^/(^P""^). 



Expansion at infinity (77-plane): Again we compute the derivative 

dr) ~ 2{p - 1) 

Tn order to use eq. (B.2) we set = 1 + /. In the new variable 

_3p+l 3p+l 2 

r] 3p-l = f 3p-l _|_ f 3p-l ^ 

and for the derivative of $ we find, using that y = {1 + /)2/(3p+i) = ^-2/(3p-i)y2/(3p-i)^ 

P 2(p-l) 2(p-l) 
= 77 3p— 1 ^ 3p— 1 

d?7 2(p - 1) ' 

3p+l 3p+l 2fT)-l") 

Applying (B.2) with /j, = rj Sf-i , A = 1, 2; = / ^p-i ^ a = and 7 = — we get 



(B.9) 



(B.IO) 



(B.ll) 



d^ 



P 



00 p 



dr) 3p + 



n=0 



2n _ 2(p-l) 
3p+l 3p+l 

p+3 3p-l 
3p+l ' 3p+l 



n 



(_l)"^r. 



(B.12) 



and finally integrating this expression with the constant $(0) = —{3p + l)/(2(p + 1)), 



00 p 



2{p+l) ' 3^,+ ! J^^ir 



2n 2p 



3p+l 3p+l 



4p+2 _ 3p-l 
3p+l 3p+l 



n 



(-l)V 



(B.13) 



3p+l 



The branches of $(77) for finite are in relation with the branch points of y{ri), solutions of y 2 = 
-i(3p-l): 



yfe = e''^^3p+T |i(3p_i)|3p+T, k = 0,±l,..., 
which gives for the positions of the branch points in the 77 plane 



(B.14) 



(B.15) 
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